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ABSTRACT 


Using  a  particular  model  of  the  photographic  point  spread  function 
derived  previously,  the  images  which  would  be  expected  from  a  reentry 
body  plus  a  luminous  wake  and  from  a  point  target  moved  uniformly  during 
exposure  are  calculated  and  compared.  The  former  images  exhibit  a 
pronounced  head  and  faint  tail  which  do  not  appear  in  the  latter. 

Two  calibration  procedures  for  images  of  the  latter  type  may  be 
considered,  based  on  the  use  of  either  the  image  growth  law  for  line 
calibration  images  or  the  image  growth  law  for  point  calibration  images. 
The  accuracy  of  either  method  depends  on  the  point  spread  function  of  the 
camera.  In  general  the  first  procedure  is  accurate  for  long  images  but 
may  have  a  50%  error  for  images  having  small  length-to-width  ratios.  The 
second  may  have  a  30%  error  for  long  images  but  is  accurate  for  short 
images;  it  is  felt  that  this  procedure  will  be  more  satisfactory  in 
practice . 
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I.  INTRODUCTION 

In  the  field  of  reentry  measurements,  photographic  images  of  line 
targets — that  is,  targets  which  are  for  practical  purposes  extended  in 
only  one  direction — are  of  considerable  importance,  since  not  only  do 
such  images  constitute  the  output  of  ballistic  cameras,  but  they  are 
produced  by  cine  cameras  when  tracking  is  imperfect  or  the  reentry  body 
has  a  wake.  When  making  reentry  measurements  with  cine  cameras,  it  is 
necessary  to  distinguish  between  the  latter  two  types  of  images,  and  to 
have  adequate  procedures  for  extracting  quantitative  data  from  each. 

In  this  memorandum  a  consideration  of  these  requirements  is  based 
on  calculations  of  the  images  which  would  in  theory  be  obtained  with 
given  camera  point  spread  functions.  The  first  part  compares  images  which 
result  from  uniform  rectilinear  smearing  with  the  images  of  reentry  bodies 
with  wakes.  The  second  part  discusses  two  calibration  procedures  for  the 
former  images.  It  will  be  recognized  that  in  such  a  treatment  many  factors 
such  as  the  Eberhardt  effect,  reciprocity  failure,  variation  of  camera 
resolution  over  the  image  plane,  and  collimator  resolution  are  neglected. 
However,  the  calculations  here  should  be  of  value  in  suggesting  further 
experimentation  and  calculation. 
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II.  IMAGES  OF  MOVING  TARGETS 

We  consider  a  photographic  system  having  a  point  spread  function 
C(x,y)  and  a  line  spread  function  A(x)  with  the  following  usual  relations 
holding: 


J  j C(x,y)  dx  dy  =  1 


j  C(x,u)  du  =  A(x) 


Consider  the  target  image  produced  by  a  moving  unresolved  point.  If  the 
power  available  from  the  point  is  W^,  and  it  follows  the  path  x  =  x(t), 
y  ■  y(t) ,  the  energy  distribution  produced  in  the  emulsion  in  exposure 
time  At  starting  at  time  0  is 


At 

^x.y)  =  W  /  C[x  -  x(t),  y  -  y(t)]  dt 


If  x(t)  =  0,  and  y(t)  =  v^t  -  and  vAt  is  set  equal  to  2a  so  the 

geometrical  image  of  the  point  moves  along  the  y  axis  from  y  =  -a  at  time 
zero  to  y  =  +a  at  time  At,  then 

At 


((x,y)  =  W  J  C(x,y  -  vt  +  a)  dt 


Changing  variables  gives 

Ep(x,y)  -  Wp  At  j  C(x,u)  du 


v+a 


(1) 


y-a 
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The  point  spread  function  of  various  cameras  can  be  approximated 
by  the  function^ 


C(x,y) 


.2 (A  -  1) 


,  2  .  2  2, 

(e  +  x  +  y  ) 


(2) 


where  A  and  e  are  characteristic  of  the  particular  camera.  The  point 
spread  function  of  the  TRAP  III  High  Speed  Camera  is  described  by  this 
expression  with  A  =  1.5  and  e  =  13  microns.  For  A  =  1.5, 


W  At  e 
E  (x,y)  =  — E. 

F 


y+a 


4air 


/ 


du 


y-a 


2  2  2 
(e  +  x  +  u  ) 


3/2 


At  the  center  of  the  image, 


E  (0,0)  =  W  At 
P  P 


1 


The  shape  of  the  image  may  be  described  by  giving  the  contour  lines 
along  which,  for  a  constant, 

E  (x,y)  =  aEp(0,0) 

-2  -3 

Contours  have  been  computed  for  ot  =  10  and  a  =  10  and  three  values 
of  a,  and  the  results  are  shown  in  Fig.  1. 
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100  MICRONS 


Figure  1. 


Equal  Intensity  Contours  for  Three  Line  Targets:  Contours 
-2  -3 

Drawn  at  10  and  10  Times  Peak 
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III .  IMAGES  EXHIBITING  RESOLVED  MAKES 

Now  consider  the  structure  of  elongated  images  when  these  are  not 
due  to  image  motion  during  exposure,  but  rather  result  from  a  small  body 
plus  a  wake.  It  will  be  supposed  that  the  geometrical  image  of  the  body 
is  a  spatial  6-function  (a  true  point)  and  that  the  geometrical  image  of 
the  wake  is  a  single  line  behind  the  body.  The  power  per  unit  length 
received  from  the  wake  will  be  supposed  to  be  a  maximum  at  the  body,  and 
to  decrease  linearly  with  distance  behind  the  body  until  it  reaches  zero. 
Thus  it  is  assumed  that  the  body  is  unresolved,  and  that  the  wake  is  not 
resolved  in  width,  but  is  resolved  in  length. 

The  energy  distribution  in  the  image  plane  of  the  optical  system 
can  be  obtained  by  adding  the  energy  distribution  due  to  the  body  to  that 
due  to  the  wake,  so  these  distributions  may  be  considered  separately. 

Also,  since  the  exposure  time,  At,  is  not  of  interest  here,  it  may  be 
assumed  that  At  *  1.  Then  if  the  powei  received  from  the  body  is  and 
the  optical  point  spread  function  is  C(x,y),  the  energy  distribution  in 
the  image  of  the  body  alone,  1^ ,  would  be  simply 

h  =  WbC(x’y) 

The  assumption  about  the  wake  given  above  requires  that  the  geometri¬ 
cal  image  of  the  wake  have  the  form 

f(x,y)  =  W^(l  +  6 (x) 


for  0  y  >_  -  d.  This  represents  a  wake  extending  a  distance  d  along  the 
negative  y  axis  from  a  body  assumed  to  be  at  the  origin;  is  the  power 
per  unit  length  at  the  head  of  the  wake.  The  total  power  from  the  v;ake 
is  thus 


00  0 

J  J  W^d  +  j)  6(x)  dy  dx 

—oo  — d 


W  ?  W  d 

r  /( d  +  y)dy  =  -f- 

-d 
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The  energy  distribution  in  the  image  cf  the  wake  alone  is 

y+d  » 

I  (x,y)  =  f  f  C(C,n)  W  (1  +  ]L-=-11)  6(x  -  0<U  dn 

u  j  j  o)  a 

n-y 

y+d 

=  f  C(x,n)  W  (1  +  dn 

J  o)  d 

n=y 


or 


Iu(x,y) 


W 

0) 

d 


d 

J  uC(x,d  +  y 
0 


u)  du 


It  will  now  be  assumed  as  before  that  the  point  spread  function  in  question 
is 


C(x,y)  -  ^ 


,  2  .  2  2, 

(e  +  x  +  y  ) 


3/2 


With  this  expression  the  image  of  the  body  alone  has  the  form 


eW 

Vx’y)  "  2 IT 


2  .  2  ^  2 
e  +  x  +  y 


3/2 


and  the  image  of  the  wake  alone  is  given  by 


W  E 

Vx’y)  =  27d 


u  du 


0  Je2  +  x2  +  (d  +  y  -  u)2] 


3/2 
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from  the  indefinite  integral 


/ 


x  dx 


r  w 


a  +  bx  +  cx 


c  V  a  +  bx  +  cx 


.  b (2cx  +  b) 
4ac  -  b^ 


we  obtain 


Xu(x,y)  - 


W  e 

U) 


2  2 

2Trd(e  +  x  ) 


2  2  2 
e  +  x*  +  (d  +  y) 


1/2  _  c2  +  x2  +  v2  +  dy 

2  2  2,1/2 
e  +  x  +  y 


The  energy  density  in  the  overall  image  is,  of  course, 


Ih(x,y)  +  I  (x.y) 

D  CO 


To  afford  an  insight  into  the  appearance  of  the  image  of  the  body 
and  wake,  a  particular  example  has  been  considered.  It  was  supposed 
that  t  =  12.5  microns  (more  convenient  for  calculation  than  the  13  microns 
determined  for  the  high  speed  camera),  that  d,  the  length  of  the 
geometrical  image  of  the  wake,  was  1000  microns,  and  that  the  total 
energy  received  from  the  body  was  ten  times  that  received  from  the  wake. 
Thus 


W  =  2  x  10-4  W, 
u  b 

and  the  energy  density  at  the  origin 

W, 

I.  (0,0)  +  I  (0,0)  = 

D  CO  0  _  2. 


is  obtained  from 
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Contour  maps  have  been  plotted  for  the  body  alone,  the  wake  alone,  and  the 
body  plus  wake  (Figs.  2,  3,  and  4).  In  each  case  the  four  contours  pass 
through  points  where  1^,  1^,  or  1^  +  1^  is  a  given  fraction  times  the 
intensity  at  the  origin,  1^(0, 0)  +  1^(0, 0) .  The  contours  can  be  considered 
to  be  the  edges  of  the  images  which  would  be  seen  as  the  power  received 
from  the  target  increased  by  multiples  of  10  from  a  level  sufficient  to 
carry  the  emulsion  well  into  saturation  at  the  center  of  the  geometrical 
image  of  the  body. 

It  will  be  noted  that  the  body  plus  wake  images  differ  significantly 
from  the  wake  images,  or  from  those  due  to  uniform  motion,  in  having  a 
definite  well-defined  head  followed  by  a  much  fainter  tail.  The  images 
shown  look  much  like  actual  resolved  images  of  large  reentry  bodies 
with  wakes. 
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IV .  CALIBRATION  AND  DATA  REDUCTION  PROCEDURES 

We  now  consider  two  techniques  for  the  reduction  of  data  from 
integrating  cameras  when  images  are  blurred  through  uniform  rectilinear 
motion  of  the  image  during  exposure.  The  methods  discussed  apply  to 
both  saturated  and  unsaturated  images,  but  are  especially  valuable  for 
the  former.  The  first  procedure  consists  of  the  following  four  steps: 

1.  Obtain  the  usual  point-source  image-growth-law  calibration 
curve  giving  image  diameter  versus  irradiance  with  a  small 
source  aperture  in  the  collimator.  (This  will  be  necessary 
in  general  for  the  measurement  of  unsmeared  target  images.) 

2.  Measure  the  width  2x^  and  the  length  2y^  of  the  elongated 
image  at  the  same  fixed  density  level  used  in  deriving  the 
above  calibration  curve. 

3.  Use  2x^  and  the  calibration  curve  to  obtain  an  initial 
approximation  to  the  target  irradiance. 

4.  Multiply  this  initial  approximation  by  the  ratio  y^/x^  to 
obtain  the  final  approximation  to  target  irradiance. 

The  second  calibration  procedure  which  has  been  investigated  is 
based  on  the  use  of  a  slit  source  aperture  in  the  collimator.  This 
procedure  is  as  follows: 

51.  Obtain  a  second  image-growth-law  calibration  curve  giving  line 
image  width  versus  irradiance  (per  unit  image  length)  with  a 
slit  source  aperture  in  the  collimator. 

52.  Measure  the  width  2x^  and  the  length  2y^  of  the  elongated 
target  image  at  the  same  fixed  density  level  used  in  deriving 
this  calibration  cu  ve. 

53.  Use  2x^  and  the  calibration  curve  to  obtain  an  approximation 
to  the  target  irradiance  per  unit  length. 

54.  Multiply  this  approximation  by  2y^  to  obtain  a  final 
approximation  to  target  irradiance. 
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A.  ANALYSIS  OF  FIRST  CALIBRATION  PROCEDURE 

Suppose  that  we  have  a  photographic  system  having  a  circularly 
symmetrical  point  spread  function  C(x,y).  Calibration  images  are  made 
with  an  unresolved  point  source,  and  the  diameters,  2rc,  of  these  images 
are  measured  at  a  fixed  threshold  level  Et>  Since  the  energy  distribution, 
Ec(x,y),  in  the  calibration  images  due  to  such  a  source  which  is  providing 
power  during  exposure  time  At  is 

E£ (x,y)  =  Wc  At  C(x,y) 

the  relationship  between  and  r ^  can  be  obtained  from 

E.  =  Wo  At  C(r  ,0) 
t  c  c 


and  hence  is 


W 

c 


At  C(r  ,0) 
c 


(3) 


This  is  the  system  image  growth  law,  which  is  thus  proportional  to  the 
reciprocal  of  the  point  spread  function. 


For  the  elongated  target  images,  x^  and  y^  can  be  calculated  from 


Eq.  (1) 


E  =  E  (x.  ,0)  =  W  At  f  C(x.  ,u)  du 

t  p  1  p  2a  J  1 


(A) 


-a 


y:+a 


Et  ■  EI,<0'yi)  •  "P  E  k  /  c(0'u)  du 


(5) 


yra 
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Also,  letting  a  be  the  ratio  between  Et  and  the  energy  density  in 
the  center  of  the  image. 


!t  ‘ » V*  s  / 


du 


(6) 


-a 


The  suggested  procedure  requires  that  2x^  and  2y^  be  measured,  a 
corresponding  be  determined  by  substituting  x^  for  rc  in  Eq.  (3),  and 

this  value  of  be  multiplied  by  y^/x^  to  obtain  a  final  computed  value 

of  target  power.  Denoting  the  computed  value  of  target  power  by  [Wp] comp 
we  have,  using  Eqs.  (3)  and  (4) 


W 


pj  comp 


At  C(xx,0) 


or 


f  C(x  ,u) 

w  i  J 

p  j  comp  vl  1_  _ 

W  x,  2a  C(x, ,0) 

pi  1 


du 


(7) 


where  y^  and  x^  are  to  be  obtained  from  Eqs.  (4)  and  (5). 

It  will  be  observed  upon  inspection  of  Eqs.  (4)  through  (7)  that 

in  addition  to  a  and  the  form  of  C(x,y),  either  Et  or  a  must  be  specified 

in  order  to  calculate  the  relative  error.  For  saturated  images  a  may  in 

-2  -5 

practice  be  between  10  and  perhaps  10  ;  the  exact  value  is  not  expected 

to  produce  a  great  effect  in  the  nature  of  the  calculated  results.  In 
this  memorandum  only  one  or  two  values  of  a  have  been  selected  for  each 
a  and  C(x,y)  in  order  to  reduce  the  computation  time  required.  With 
this  limitation,  the  accuracy  of  the  proposed  calibration  method  has  been 
evaluated  for  three  point  spread  functions  having  the  form  given  by 
Eq.  (2),  with  values  of  1.5,  2,  and  2.5  for  A;  these  cover  the  range  of 
image  growth  laws  encountered  so  far.  The  calculated  values  of  the 
relative  accuracy  are  shown  in  Fig.  5  for  various  choices  of  the  parameters. 
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0  5  10 

IMAGE  LENGTH 
IMAGE  WIDTH 

Figure  5.  Ratio  of  Computed  Irradiance  to  Actual  Irradiance  as  a  Function 
of  the  Length-to-Width  Ratio  of  Smeared  Images  for  First  Cali¬ 
bration  Procedure.  (A  is  One-Half  the  Limiting  Slope  of  the 
Image  Growth  Law  on  a  Log-Log  Plot.) 

It  can  be  seen  that  the  proposed  procedure  is  fairly  accurate  for 
images  which  are  only  slightly  elongated,  but  may  introduce  an  error  as 
grea*-  as  33%  for  very  long  images.  This  is  considered  satisfactory,  since 
in  practice  the  images  to  be  measured  are  generally  not  very  long,  and 
errors  of  this  magnitude  are  usually  acceptable. 

B.  ANALYSIS  OF  SECOND  CALIBRATION  PROCEDURE 

The  alternate  calibration  method  considered  for  such  targets  is 
based  on  using  a  slit  source  in  the  collimator.  If  the  power  per  unit 
length  in  the  geometrical  image  of  the  slit  is  Wg,  the  exposure  time  is 
At,  and  the  geometrical  image  of  the  slit  coincides  with  the  y  axis,  the 
energy  distribution  per  unit  length  in  the  calibration  image  is 

CO 

E  (x,y)  =  W  At  f  C(x,y)  dy  =  W  At  A(x) 
s  S  J  s 
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where  A(x)  is  the  system  line  spread  function.  If  a  series  of  such  images 
is  made  with  different  values  of  Wg  and  the  half-width  xg  of  each  image 
is  measured  at  a  fixed  level  Et>  the  resulting  image  growth  law  is 
obtained  from 


E 


t 


W  At 
s 


j  C(xs,y)  dy 


and  hence  is 


W 

s 


At  j  C(xg,y)  dy 


(8) 


In  order  to  determine  responsible  for  an  image  of  the  form  given  by 
Eq.  (1),  some  estimate  of  the  effective  image  length  must  be  obtained 


from  the  image. 

of  W  ,  Iw 
P  I  T 


p| comp  2* 


If  this  estimate  is  denoted  by  2  a,  the  computed  value 
is 


W 


pi  comp  2 


2  a  W 


Hence,  using  Eqs.  (4)  and  (8), 


W 

p  comp 
W 

P 


2 


2  a  W 

s 


W 

P 


I  /  C(x1>y)  dy 

-a 


j  C(x1,y)  dy 


%  % 

The  only  reasonable  choice  for  a  appears  to  be  a  =  where  y^  satisfies 
Eq.  (5).  The  choice  a  =  y^  -  x^  which  might  be  considered  leads  to  the 
unsatisfactory  result 

[w  I 

Lim  |  p| comp  2  _  n 
a-0  W 

P 
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With  a  =  we  have 

a 

„  |  hf  c(v’r)  d* 

P)  com£.  2  =  -a 

W  y  1  - 

p  1 

J  CCx^y)  dy 


We  note  that  as  a  increases 

Lim  [Wplcomp  2  _ 
a-*°°  W 

P 

and  as  a  decreases 

T  .  I W  1  „  2y.  C(x. ,0) 

Lim  1  p|  comp  2  _  J 1 _ 1 _ 

CO 

J  C(x1,y)  dy 


a-*-0  W 

P 


(9) 


Comparing  Eq.  (9)  with  Eq.  (7),  we  obtain 


|wl  2 

l  pi  comp  2 

|w  1 

l  pi  comp 

’  ' 

2x.  Cvx  ,0) 

W 

W 

OO 

P 

P 

(  C(xx,y)  dy 

— oo 


A  _  l  c2(A-1) 

When  C(x,y)  =  — -  — r - - - r — -  it  is  found  that 

7T  .1,1,  A 

(e  +  x  +  y  ) 


f  C (x,y)  dy 


r(i/2)r(A  -  i/2) 
r(A) 


A  -  1 


2 (A-l) 


2  2 
e  +  x 


A-l/2 
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and  hence 


w 

p  comp  2 

W 

p  comp 

2T  (A) 

xi 

w 

p 

w 

p 

r(i/2)r(A  -  1/2) 

2  2  1/2 
e  +  x^ 

This  enables  the  relative  error  for  calibration  with  a  slit  source  to  be 
obtained  from  the  calculations  which  gave  Fig.  5;  the  results  for  the 
three  values  of  A  considered  are  shown  in  Fig.  6. 


Again,  the  accuracy  of  this  procedure  depends  on  the  system  image 
growth  law.  It  will  be  noted  that  this  procedure  gives  the  correct 
irradiance  for  large  values  of  the  ratio  y^/x^,  t>ut  not  as  accurate  when 

y^/x^  is  close  to  unity.  This  procedure  is  considered  inferior  to  the 
previous  one  because  most  elongated  images  are  relatively  short  and  thus 


0  5  10 

IMAGE  LENGTH 
IMAGE  WIDTH 


Figure  6.  Ratio  of  Computed  Irradiance  to  Actual  Irradiance  as  a 
Function  of  the  Length-to-Width  Ratio  of  Smeared  Images 
for  Calibration  with  Slit.  (A  is  One-Half  the  Limiting 
Slope  of  the  Image  Growth  Law  on  a  Log-Log  Plot.) 
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fall  in  the  region  where  the  previous  method  is  more  accurate,  and 
because  the  second  method  requires  an  additional  set  of  calibration 
measurements  while  the  first  makes  use  of  the  standard  calibration  curve 
needed  for  data  reduction  of  unsmeared  targets. 
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